Quantum metrology at the Heisenberg limit with ion traps 



o : 
o . 

(N ■ 

W> 
=5 

<■ 

On i 



> 

(N . 
00 ■ 

o ■ 
oo : 
o . 
\o ■ 
o : 

^ ■ 

Q-i! 
^— > ■ 
C ■ 

=5 

cr 



X 



D.A.R. Dalvit, 1 R. L. de Matos Filho, 2 and F. Toscano 3 ' 2 

1 Theoretical Division, MS B213, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
z Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68.528, 21.941-972, Rio de Janeiro, Brazil 
s Fundagao Centro de Ciencias e Educagao Superior d Distancia do 
Estado do Rio de Janeiro, 20943-001 Rio de Janeiro, RJ, Brazil 
(Dated: February 1, 2008) 

Sub-Planck phase-space structures in the Wigner function of the motional degree of freedom of a 
trapped ion can be used to perform weak force measurements with Heisenberg-limited sensitivity. 
We propose methods to engineer the Hamiltonian of the trapped ion to generate states with such 
small scale structures, and we show how to use them in quantum metrology applications. 
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I. INTRODUCTION 



The determination of small parameters has recently ac- 
quired a substantial improvement through quantum mea- 
surement, as it is now known that using probes prepared 
in judiciously chosen quantum states can increase their 
sensitivity to perturbations. As a consequence, quantum 
metrology has become a subject of great practical inter- 
est 0. The estimation of an unknown parameter of a 
quantum system typically involves a three-step process: 
the initial preparation of a probe in a known quantum 
state, the interaction between the probe and the sys- 
tem to be measured, and a final read-out stage, where 
the state of the probe is determined. Typical situations 
are those when the system imprints an unknown param- 
eter x onto the probe through a unitary perturbation 
U x = cxp(-ixG), where the generator G is a known 
Hermitian operator. The unknown small parameter x of 
the system can be inferred by comparing the input and 
output states of the probe. This framework captures two 
important tasks in quantum metrology: i) high precision 
phase measurements x = 9, where Ug generates a rota- 
tion in phase space of the quantum state of the probe 
around the origin, and ii) detection of weak forces that 
induce a linear displacement x = s of the quantum state 
of the probe in some direction of phase space. 

The accuracy of the parameter estimation is limited 
by the physical resources involved in the measurement. 
Techniques involving probes prepared in quasiclassical 
states, such as coherent states of light, have sensitivities 
at the standard quantum limit (SQL), also known as the 
shot-noise limit in the phase detection situations. Indeed, 
in the usual dimensionless phase-space used in quantum 
optics the Wigner distribution of a coherent state, with 
n mean number of photons, is centered at a distance 
~ \ffi from the origin, with a width 1/2. Thus, the 
associated input and output states are distinguishable 
(approximately orthogonal) when their respective phase- 
space distributions are displaced by a minimal distance 
1/2 ~ O(n ) from one another, i.e., the SQL for weak dis- 
placement measurement does not depend on the number 



n of photons involved in the measurement. For phase 
detection the smallest noticeable rotation occurs when 
the centers of the phase-space distributions of the input 
and output states have an angular distance (1/2)/a/H 
measured from the origin, i.e., the SQL for phase mea- 
surement scales as l / \ffi. 

Using the same physical resources in addition to quan- 
tum effects, such as entanglement or squeezing, sub SQL 
precision can be achieved @, IE & IE IE IE IE • For exam- 
ple, probes prepared in quantum superpositions of, say, 
coherent states of light, have Wigner distributions with 
sub-Planck phase-space structures of typical linear size 
IE- ^ was recently shown [ljj how these struc- 
tures can be used to achieve Heisenberg-limited sensitiv- 
ity in weak force and phase measurements: the approx- 
imate orthogonality that allows to distinguish the input 
and output quantum states of the probe occurs when the 
peaks (valleys) of the sub-Planck structures of one come 
on top of the valleys (peaks) of the sub-Planck structures 
of the other. The minimal linear displacement required 
for this destructive interference is ~ 1/V^, defining the 
so-called Heisenberg limit (HL) for weak force detection. 
For rotations the sub-Planck phase-space structures of 
the Wigner function of the input state have to be at a 
typical distance ~ from the origin, thus the minimal 
rotation angle is ~ 1/n, that defines the HL scale for 
phase detection. Hence, the linear size of the sub-Planck 
structures sets the sensitivity limit on a probe and states 
that saturate the limit on the smallest phase-space struc- 
tures can allow one to attain Heisenberg-limited sensi- 
tivity. Sub SQL sensitivities, approaching the ultimate 
Heisenberg limit, have been recently achieved experimen- 
tally using multi-qubit states formed by entangled inter- 
nal de gree s of freedom of photons (polarization entangle- 
ment) ^lUElEEEil] , an d °f trapped ions (spin entan- 
glement) |l5l llo . Il7j . 

Sub-SQL precision measurements can also be achieved 
with continuous variables, such as the amplitude/phase 
of a photon trapped in a QED cavity, or the center of 
mass motion of trapped ions [Ts[ . Such spatial modes, 
that can be approximately described as harmonic oscil- 
lators, can be prepared in nonclassical quantum states, 
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e.g. superpositions of M coherent states, that possess 
sub-Planck phase-space structures. In |ldj | it was shown 
how to use them for Heisenbcrg-limitcd precision mea- 
surements of weak forces and phase measurement, by en- 
tangling the oscillator with a two-level system. It was 
also proposed a concrete implementation of this strat- 
egy in the context of cavity QED (that can be extrap- 
olated to trapped ions): a two- level atom gets entan- 
gled with the cavity mode through the (resonant or non- 
resonant) Jaynes-Cummings interaction in such a way 
that, when the atom reaches the center of the cavity, the 
photon field is a quantum superposition of two coherent 
states (M = 2). Such "cat state" has sub-Planck oscilla- 
tions parallel to the line joining the two coherent states, 
and is therefore Heisenberg-limited sensitive to perturba- 
tions that induce displacements that are perpendicular to 
that line. After the perturbation is applied, the Jayncs- 
Cummings evolution is undone, and the populations of 
the two-level atom are measured once it leaves the cavity. 
This provides information about the perturbation, that 
can be estimated with HL precision. An analysis of the 
decoherence processes ^} that may affect this scheme 
shows that this proposal should be within reach of cavity 
QED and ion trap experiments. 

One disadvantage of the M = 2 states is that their 
sensitivity gradually degrades as the direction of the per- 
turbing force moves away from the direction orthogonal 
to the line joining the two coherent states. Higher or- 
der (M > 2) superpositions of coherent states on a circle 
do not suffer from these limitations. In the context of 
cavity QED, there have been a number of schemes pro- 
posed to generate such states. In some of these proposals 
jTol I2H EH I22I it is necessary to send a sequence 
of atoms through the cavity and perform conditional 
measurements (post-selection), limiting the probability 
of succeeding in engineering the cavity in the desired co- 
herent state superpositions. There are also schemes that 
involve only one atom, and make use of a single-atom 
interferometric method p3 |. or a dispersive interaction 
between the atom and the cavity field [25|. The first 
scheme has the disadvantage of requiring post-selection, 
while the second one needs very long interaction times. 
Similar strategies have been proposed in the literature in 
the context of trapped ions |2jj, |2j| ■ 

In this paper we propose ways for deterministic gen- 
eration of higher order (M > 2) coherent state superpo- 
sitions with trapped ions, and how to use them to mea- 
sure weak forces at the Hciscnbcrg limit. These states 
do not suffer from the limitations of the M = 2 states 
generated via the linear Jaynes-Cummings interactions, 
as discussed above. In particular, we will concentrate on 
how to prepare probes with M = 4 ("compass state") 
and how to use them to measure small displacements 
(weak force detection). Our setups can also be used to 
measure small rotations in phase space (phase detection) 
by simply adding an appropriate displacement before the 
application of the perturbation, as described in [lflj . 

The paper is organized as follows. In Section II we 



discuss the sensitivity to perturbations of the M = 4 
compass state and we review the strategy for weak force 
measurements. In Section III we review the engineering 
of the ion-laser interaction for a single trapped ion. Sec- 
tion IV contains two possible approaches for generating 
compass states: the first one uses conditional quantum 
gates between the internal and the motional degrees of 
freedom of the ion, and the second one uses a nonlin- 
ear ion-laser engineered interaction. In Section V we 
implement the scheme for measuring weak forces with 
Heisenberg-limited sensitivity. Finally, Section VI sum- 
marizes our central results. 



II. HEISENBERG-LIMITED QUANTUM 
METROLOGY WITH CONTINUOUS 
VARIABLES 

In this section we review how to perform continu- 
ous variable based quantum metrology at the Hciscnbcrg 
limit by using sub-Planck phase-space structures present 
in circular coherent states of a harmonic oscillator (the 
probe) [13 ■ In the following wc will focus on perturba- 
tions that induce a small phase-space linear displacement 
of the quantum state of the probe, for example when a 
small force is exerted onto the probe. The unitary oper- 
ator describing such a displacement is D{(3) = e 13 ^^ 13 a , 
where (3 is an arbitrary small displacement in phase-space 
with a magnitude <C \a\. In this approximation the 
unitary operation D(($) takes any coherent state \a) to 
e iim(/3a*)| a + e 2Um ^ a ^\a). Perturbations inducing 

small rotations R(0) = e l9aia (9 <C 1) in phase space can 
be treated in a similar way, as described in [Toj . 

Circular coherent states are a special kind of states 
of harmonic oscillators, formed by the superposition of 
M coherent states equidistantly placed on a circle C of 
radius |a|. They are defined as 

M M 

|cat M ) = 7 =^e^ \e^a), (1) 

V M k=l 

where ip^ = 2nk/M, 7^ are arbitrary phases, and M = 
l + C(e _ l Q l ) is a normalization constant. When \a\ > 3, 
it is a good approximation to take Af = 1, and we will 
assume this is the case in the following considerations. 
These non-classical states of a harmonic oscillator have a 
mean number of excitations n = (catju |a^a|catA/) m |a| 2 . 
Examples of these states are the Schrodinger cat state 
(M = 2), considered in detail in [T^j f° r weak force and 
phase measurements, and the compass state 

|cat 4 ) = i (e 171 \ia) + e* 72 1 - a) + e* 73 \-ia) + e* 74 \a)) , 

(2) 

that we consider here. 

Figure depicts the Wigner function of the compass 
state and shows the sub-Planck phase-space structures, 
that oscillate with a typical wavelength ~ These 
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FIG. 1: (Color) The Wigner functions in the a plane for: (a) the compass state |cat4.) with a = 3, and (b) the displaced 
compass state D{j3)\ca±4) — |cat4(s D )) where /3 — e l{p s a/\a\ is a displacement of magnitude \/3\ — s — ir/2b+(tp)\a\ in the 
direction relative to a given by tp = 7r/3. Here 6+(<p) = cos{ip) + sin(ip). The black arrow in (b) indicates the direction of the 
displacement (3. In (c) we display the product of the unpertubed and the perturbed Wigner functions. When performing the 
integration over the a plane of this product of Wigner functions, that is equal to the overlap |{cat4|cat4(s )}| 2 , the negative 
contributions (in blue) cancel the positive ones (in red), leading to quasiorthogonality. 



structures are responsible for the Hciscnbcrg-limited sen- 
sitivity of these states for quantum metrology applica- 
tions. Indeed, the smallest linear displacement needed to 
distinguish the unperturbed |cat4) from the perturbed 
|cat4(s)) = D{j3)\cati) state of the probe (and, there- 
fore, to attain quasi-orthogonality) is s ~ l/|q|, for any 



direction (3 = e lv sa/\a\. As it was shown in 



the sub- 



Planck structures of the Wigner functions of the circular 
states determine the oscillatory behavior of the fidelity 
function f(s) = |(cat4|cat4(s))| 2 , that for the compass 
state reads 



/(») * i 



fe=l 



= e" s cos 2 (b + (tp)\a\s) cos 2 (b_((p)\a\s) , (3) 

where ctk = e lVk a and b±(<p) = cos(y) ± sin(</?) (cp is 
the angle between (3 and a). Here we have neglected 
contributions (ai\ak + (3) ~ C(e~' Q ' ) for I ^ k, that is 
a good approximation for \a\ > 3, and we have used that 
(a k \a k + (3) = e~^ 2 exp{ilm(f3a* k )}. When e^ 2 w 1, 
valid for |/3| = s <C 1, we can re- write the fidelity function 
as 



(b-(<p)\a\s). (4) 



Note that /(s) does not depend on the phases jk that 
enter in the definition of the compass state |cat4) in 
Eq.(j2J. We can see in Figure the oscillatory behav- 
ior of the fidelity f(s) as a function of the magnitude 
of the linear displacement s = \/3\ with a typical wave- 
length ~ l/M- The minimal displacement to achieve 
quasi-orthogonality, s , can be obtained equating to zero 
the approximation in Eq.QJ because we see in Fig[21 



that this approximation well describes the fidelity beyond 
s . Thus, the minimal displacement to achieve quasi- 
orthogonality is approximately given by 



(5) 



2%)M 



where b(ip) = max{|& + (<p)|, | b-(tp) \} (sec Fig.©). The 
scaling of s D as l/|a| implies that one can detect dis- 
placement perturbations with Heisenberg-limited sensi- 
tivity. Note that b+(ip) and b-((p) are never simultane- 
ously zero in the interval [0, 2tt], which means that for 
any direction of the displacement f3 there is a minimum 
(finite) displacement for which the unperturbed and the 
perturbed compass states are quasi-orthogonal. Contrary 
to the cat state (M = 2), whose HL sensitivity degrades 
as the direction of the perturbation /3 tends to the di- 
rection of a (i.e., as — > 0), compass states arc HL 
sensitive for all directions if of the perturbing force. In 
a similar way, one can show that compass states have 
Heisenberg-limited sensitivity to rotation perturbations, 
that is, the minimal rotation angle that can be detected 
is6*~s/|a|~l/|a| 2 . To achieve this sensitivity for ro- 
tation perturbations, the circular coherent state has to 
be first translated in phase space so that the displaced 
circle C contains the origin of phase space |Hj . 

Once the direction of the perturbation (3 is known, i.e., 
the angle (p), the information about its magnitude is en- 
coded in the fidelity function /(s) between the unper- 
turbed and the perturbed states of the probe. In order 
to measure the fidelity, a measurement strategy was pro- 
posed in ^(j j that is summarized in Fig. It involves 
entangling the probe with an ancillary system, a two-level 
system (TLS), and designing their interaction U in such 
a way that the information about the overlap between 
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FIG. 3: Circuit diagram to measure a small perturbation U x 
at the Heisenberg limit. The probe is a harmonic oscillator 
and the ancilla a two-level system. The entangling evolution 
U is such that |(*1/| < I > )| 2 ~ |(cat4|cat4(a:))| 2 , where |cat4) is the 
compass state defined in Eq.^, and |cat4(a;)) = U x \ca.t4) the 
perturbed compass state (see text for details). 



FIG. 2: The fidelity f(s) = |(cat4|cat4(s))| 2 as a function of 
the magnitude of the linear displacement s = \/3\ (full line) 
where /3 = e llp sot/\a\ with <p = 7r/3. The circles (o) corre- 
spond to points given by the approximation in Eq.©. The 
dashed line is the approximation in Eq.J^J The displacement 
corresponding to quasiorthogonality, i.e., s in Eq.lJKJ is indi- 
cated in the plot. 



the unperturbed and perturbed probes states can be in- 
ferred from measurements of the populations of the TLS. 
Starting with an initial joint state \^>i) = |0) <X> \g), where 
|0) is the vacuum state of the probe and \g) the lower 
state of the ancilla, the evolved states, with and with- 
out the perturbation applied, are |$) = U x U\^>i) and 
\^>) = U\^i) respectively. Here U x denotes the unitary 
operator corresponding to the perturbation. The inter- 
action evolution U has to be designed in such a way that 
|(>3>|$)| 2 ps |(cat 4 |cat 4 (.T))| 2 . Then, we have to undo the 
unitary evolution U to obtain a final entangled state of 
the composite system 



(6) 



where P e and P g = 1 — P e are the probabilities of mea- 
suring the TLS in levels e and g, respectively, given by 



\(^i\^f)\ 2 |(cat A/ |cat A/ (a;))p 
l(0|*!>| 2 * 



l<0|*|)| 2 



(7) 



As mentioned in [lj}, this strategy can also be used to 
measure the Loschmidt echo, which quantifies the sensi- 
tivity of a quantum system to perturbations. Also, sim- 
ilar hybrid continuous variable-qubit systems have been 
proposed for quantum information processing and quan- 
tum computation [28L 1291 . 

Later in the paper we will show how to create com- 
pass states of the vibrational degree of freedom of a 
trapped ion, and how to measure weak perturbations at 
the Heisenberg limit by entangling the motional degree 
of freedom with a two-level system, corresponding to two 
internal hypcrfine levels of the ion. 



III. ENGINEERING THE ION -LASER 
INTERACTION IN AN ION TRAP. 

Let us consider a single ion confined inside a harmonic 
trap. In good approximation the quantized motion of 
the ion center-of-mass along each spatial dimension can 
be described by a quantum harmonic oscillator. When 
the ion is illuminated by laser light quasi-resonant to 
one of its electronic transitions its motional degrees-of- 
freedom can be coupled to the electronic ones via photon- 
momentum exchange. The laser excitation can be done 
in several different ways, giving rise to a large number of 
possible interaction Hamiltonians between electronic and 
motional degrees-of-freedom. Here we will be interested 
in two basic types of laser excitation in a situation where 
the motional sidebands can be spectroscopically well re- 
solved. Moreover, we will consider the excitation of only 
one motional degree-of-freedom. 

The first excitation scheme of interest is the Raman 
excitation of a dipolc-forbidden electronic transition be- 
tween two hypcrfine electronic states, on resonance to a 
given motional sideband. This can be done via the off- 
resonant excitation of an intermediary electronic level by 
two laser fields of adequate frequencies (see Fig.QJa)). In 
this case, the interaction Hamiltonian, in the interaction 
picture, is given by [30| : 

Hi = pinole** a+/jfc(n,»7) a k +H.c. , (8) 

where the operator- valued function fk{n,r)) depends on 
the motional number operator n = a) a, 



fk(n,v) = e 



- 2/2 E 

1=0 



(irf) 



2l+k 



l\(l + k)\ (h-l)\ 



(9) 



and we define n!/(n — l)\ = [h — (I — l)][n — (I — 2)] ... n. 
The operators &+ and a are the electronic flip operator 
between the states \g) and |e) and the annihilation opera- 
tor of a motional quantum, respectively, fig = \D,q \ e 1 ^ is 
the effective Raman Rabi frequency and k corresponds to 
the excitation of the fcth lower motional sideband (in the 
example of Fig. Ufa) k — 1). The Lamb-Dicke parameter 
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excitation (k 
as: 



0), the function fo(n, rf) can be written 



fo(n,v) ~ Ao + Axh, 



(12) 



with Aq = 1 — r/ 2 /2 and Ai = — t; 2 . Under this condition, 
the choice of (j> = in the Hamiltonian (|SJ) leads to the 
evolution operator 



|K> 



- n+1 

- n 
n-1 



-n+1 
- n 
n-1 



Rc(0) 



e -iH lt /H = e . 



(13) 



FIG. 4: Electronic (internal) and motional (external) energy 
levels (not to scale) of a trapped ion, coupled by two Raman 
lasers, (a) corresponds to the first red sideband transition 
between two electronic levels, and (b) corresponds to the first 
blue sideband transition via the virtual excitation of an elec- 
tronic level. 



77 can be defined as 77 = (<5kL • u)Axq, where Axq is the 
spread of the motional wave function in the ground state 
of the harmonic potential and u is a unit vector giving 
the direction of the motional degree-of-freedom we are 
considering. The vector <5kL = kj^ kL 2 is the difference 
of the wave vectors of the two Raman lasers. The value 
of the Lamb-Dicke parameter r\ can be changed via the 
geometry of the laser beams in two different ways. The 
first one consists in changing the modulus of the vector 
<5k. This modulus reaches its minimum value when the 
laser beams are co-propagating and its maximum value 
when the laser beams are counter-propagating. The sec- 
ond way consists in changing the direction of the vector 
(5k with respect to the unit vector u. So, via the laser 
geometry, it is possible to change the value of 77 in some 
extent. In particular, using co-propagating laser beams, 
one can reach the Lamb-Dicke regime where 77 <C 1. 

Note that when ?/ <C 1 the function fk(n, 77) « {iii) k ■ In 
this case, the excitation of the carrier resonance (k = 0) 
leads to the interaction Hamiltonian: 



Hi 



1 



h\n \e l,f ' a+ 



H.c. 



(10) 



This Hamiltonian allows one to perform rotations of the 
electronic state of the ion according to the evolution op- 
erator 



U e {v) 



-iHit/h 



- ) (v x a* 



■JyUyj 



(11) 

where 9 = |flo|i, v x = cos(</>), v y = — sin(0) and a x , 
a y are the Pauli operators. Our weak force measurement 
scheme involves reversing every unitary evolution applied 
in the process. Hence, the inverse evolution Ug(—v) can 
be implemented simply by changing by (f> — > <f> + n the 
phase of the effective Raman Rabi frequency Oo . 

In situations where the value of the parameter 7/ is not 
extremely small (77 ss (0.1,0.2)) terms of order i] 2 should 
be taken into account in Eq. @. In this case, for carrier 



where v = -|f2 |L4 /2 and 6 = — |O |t^4i/2. The action 
of this operator corresponds to a rotation in phase-space 
of the motional state of the ion conditioned to its elec- 
tronic state: 

i? c (0)|a>| T*> = e iv \ad°)\U) 
Rc{6)\a)\[ x ) = e~ iv \ae- i9 )\U), (14) 

where \a) represents a coherent motional state, | 1 X ) = 
(|e) + \g))/y/2 and | i x ) = (|e) - \g))/V2 are eigenstates 
of the a x Pauli operator, and |e) and \g) represent, re- 
spectively, the upper and lower electronic states of the 
transition driven by the lasers (eigenstates of the a z Pauli 
oparator). The inverse operation R c {—9) can be imple- 
mented choosing <f> = tt. 

The dependence of the Hamiltonian (jSJ) on the mo- 
tional number operator n is determined by the function 
fk (71,77). Some control of this dependence is obtained 
by changing the value of the Lamb-Dicke parameter 77. 
However, this control is rather limited and it is of inter- 
est to find ways of extending the possibilities of tailoring 
the n-dependence of the Hamiltonian (JSJ) . As has been 
shown in |3l| , some extra tailoring can be done by excit- 
ing the same vibrational sideband of the ion by N pairs 
of Raman lasers of arbitrary effective Rabi frequencies Qj 
and Lamb-Dicke parameters rjj (J = 1, . . . , N). Remem- 
ber that the parameter 77-,- , corresponding to a given pair 
of Raman lasers, can be controlled by varying the geom- 
etry of those laser beams. Under such excitation scheme 
the resulting interaction Hamiltonian maintains the gen- 
eral form of the Hamiltonian © . The combined effect of 
the N laser pairs appears in the operator function fk(n), 
which is transformed in the engineered function ff. (n). 
This function can be written in the form of a Taylor se- 



with the coefficients A p given by 



(15) 



-4,, 



fe-'.'/^M if p = 0, 



3=1 

N 



(-!)*>£ e-^^L^f^. if p ^o, 



(16) 
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where f2^ is a reference Rabi frequency and the coeffi- 
cients a P j are definded by: 

oo 2m 

a P ,j = E (-1)—^) ^ . (17) 

rn—p v 7 

The S'm'' are the Stirling numbers of first kind [U that 
count the number of permutations of m elements with 
p disjoint cycles. The above expression shows that, for 
a given set of parameters r/j, N of the coefficients A p 
can be independently fixed by the values of the TV Rabi 
frequencies Qj. The values of the Rabi frequencies arc 
the solutions of the set of N linear equations obtained 
from Ea. (|16|) by fixing the values for the N coefficients 
A p . The free parameters r\j can be used to optimize 
the coupling between the electronic and the vibrational 
degrees of freedom. This excitation scheme opens the 
possibility for engineering a large number of interaction 
Hamiltonians. 

In the next section we will be interested in engineering 
the function 

j fin,?]) = A + A in + A 2 n 2 + 0(r,f nax n 4 ), (18) 

where only the values of A 2 and A3 = have to be fixed 
independently. This can be done with N = 2 pairs of 
laser beams in Raman configuration. If the two laser 
pairs excite the carrier resonance (k = 0), the resulting 
interaction Hamiltonian will be 

Hi = \w L \e«&JF%Ti)+Tl.c. . (19) 

For (j> = in the above equation one gets the following 
Kerr-type evolution operator 

V(fa) = e~ i ^ Ilt / h = e - ^ * 1 e ~ i 't>^^ fl (20) 

with fa = \n L \Ajt/2 (J = 0, 1, 2). We index V just with 
fa because, when applying this operation to a trapped 
ion, only the value taken by fa will be of importance. 
Note that the evolution generated by V(fa) may be in- 
verted by setting = ir in Eq. l |T5j l . 

The second basic type of laser excitation scheme we 
are interested in is the Raman excitation of one motional 
sideband via the virtual excitation of a given electronic 
transition (see Fig. P^b)). In this excitation scheme only 
the motional degree-of-freedom is excited conditioned on 
the occupation of a specific electronic level. In the exam- 
ple of Fig. J2fb)), if the ion is in the electronic state \g) 
nothing happens, whereas if the ion is in the electronic 
state |e) its vibrational motion will be excited. In this 
case the interaction Hamiltonian, describing the action 
of the lasers on the motional degree-of-freedom, is given 

by mm 

Hi = ^H\n \ e^f k (h, 77) a k + H.c. . (21) 



When the first motional sideband is excited (k = 1) in a 
situation where ri <C 1, the above Hamiltonian simplifies 
to 

Hi = -7i|r2o|e # a + H.c. . (22) 

The time-evolution generated by the Hamiltonian (|22|) is 
equivalent to action of the displacement operator D(a) 

U{t) = e- i6lt l % = D (-\r,Po\ e~ 1 ^ . (23) 

For this reason, the Hamiltonian (|22(l can be used to 
coherently displace the motional state of the ion in phase- 
space conditioned on the occupation of a given electronic 
level. In order to stress the dependence of the action 
of the above Hamiltonian on the electronic state we will 
represent the evolution operator (|23|l by D c (a). 

IV. DYNAMICAL GENERATION OF THE 
COMPASS STATE IN ION TRAPS 

According to the general procedure to measure small 
perturbations, summarized in Fig.(PJ), the measurement 
of weak forces that couple to the motional degree of free- 
dom of a trapped ion involves the dynamical generation 
of an intermediate maximally entangled state between 
the electronic (the ancilla) and the motional (the probe) 
degree of freedom of the ion, i.e., 

I*} = ^ |cat 4 ) I T) + ^ |cati> I |) , (24) 

where | T)> I !} are two orthogonal internal states of the 
ion, and |cat4), |cat4) are compass states (Eq.@) of the 
center-of-mass motion of the ion. In the following we 
describe two ways to dynamically generate these states 
starting from an initial state of the form |\E^) = |0) ® 
\g), where |0) is the ground state of the center-of-mass 
motion of the ion and \g) is the lower electronic state 
considered. The first approach combines quantum gates 
over the internal states and conditional linear operations 
over the center-of-mass motion, described in Section III. 
The second approach uses the engineering of a non-linear 
ion-laser interaction of the Kerr-type, also described in 
Section III. 



A. First approach 

In this approach the set of unitary operations that 
leads to the state in Ea. H24|) is described in Fig. [S] The 
first type of unitary operations in the sequence are the 
quantum gates Ug{v), given by the carrier 0-pulses in 
Ea. (|ll|) that rotate the internal states of the ion around 
the Bloch vector v by an angle 9. We only perform ro- 
tations around the Bloch vector v — (0,v y ,0) that can 



7 



U 



|a = 0>- 



D c (2 ae'T) 




R(f) 




D(-a) 




D c (-i'a ) 









IMVy) 



IMVy) 



U»(v v ) 



IMVy) 



u 



a = 0>- 

[*> 











D(a) 




V(f) 












1 






k 



FIG. 5: Circuit diagram to generate the compass state via 
quantum gates. Starting with an initial state |*;) = |0) ® \g) 
of a trapped ion, where |0) is the ground state of its center- 
of-mass motion, and \g) is the lowest of the two hyperfine 
electronic states considered, this sequence of unitary opera- 
tions generates an output state of the form in Eg. 1241 . 
with | 1} = \e) and | [) = \g), and |cat4) and |cat4) compass 
states. The unitary evolution U in Fig- (O is composed of this 
sequence of operations. 



be carried out by choosing the phase of the effective Ra- 
man Rabi frequency Hq equal <p — — 7f/2. Setting the 
time duration of the pulse 9 = tt/2 = \^la\t^/ 2 we can 
implement U v /2{vy), that performs 7r/ 2- rot at ions, and 

with 6 = it = |f2o|*7r we implement U„/2(v v ), that per- 
forms 7r-rotations. The other type of unitary operations 
in the sequence of FigEl affect the motional state of the 
ion conditioned on the state of the electronic degree of 
freedom. D c (a) given in Ea. (|23[) displaces a coherent 
state of the vibrational motion of the ion if the internal 
electronic part is in the upper state |e), and does noth- 
ing if the electronic part is in the ground state. That 
is, L> c (a)|a)|e) = e tIm ( aa *)|a + a)|e), and D c (a)\a)\g) = 
\a)\g). The conditioned operation R c (n/4) in Fig. [3]per- 
forms phase space rotations of a coherent state of the 
vibrational motion of the ion according to Eq. I|14|) . 

The final entangled state \^f) of the sequence of unitary 
operations of Fig. [3] is of the form given in Eq. l|24|l . with 
| |) = |e) and | J,) = \g) . and the vibrational motion of 
the ion in compass states of the form 



cat 4 ) 



e w \ia) + e 



|cat 4 ) = ^{e w \ia) 



a 



ia)+e- w \a)) 
ia)-e~ iv \a)) 



that are equally oriented in phase space. 

A realistic estimation for the total time needed to per- 
form the unitary operations in this approach can be 
done using the experimental parameters in [32|. For 
the U^/2(v y ) operation the Raman Rabi frequency is 
of the order f2o/27r m 250 kHz, so the 7r-pulse dura- 
tion is about 1/is. This implies a total time for single 
qubit operations of about 2.5/^s. The conditional dis- 
placement D c (a) of amplitude \a\ — i]CIqt/2 can be im- 
plemented with a Lamb-Dicke parameter 77 sa 0.15 and 



FIG. 6: Circuit diagram to generate the compass state via 
engineered Kerr nonlinearity V . The input state is the same 
as in Fig. The output state <&) has the form in Eg. 1241 . 
with I f) = I 1 X ) and J.) = J. x ), and |cat4) and |cat4) com- 
pass states. The unitary evolution U in Fig.0 is composed 
of this sequence of operations. 



an effective strength of the Raman laser configuration 
flo/2ir w 300kHz. Let us suppose that we create compass 
states with amplitude \a\ ~ 3, thus for the maximum dis- 
placement in the process a = 2a we have a duration time 
r f=a 21/j,s. Therefore, the total time for the conditional 
displacements in the process is aproximately 42/is. For 
the duration of the conditional rotation R c {9) we have 
6 = it/4 = iloi^t^ 1 4/2, and using also rj m 0.15 and 
Qq/2tt sa 300kHz, we get t„./ 4 « 27/ns. Hence, the total 
time needed to generate the state |\&) is approximately 
72/^s. 



B. Second approach 

In this approach the set of unitary operations that 
leads to the state in Ea.(|24|> is described in Fig. 
The displacement operator D c (a) acting on the vibra- 
tional state of the ion can be implemented along the 
lines of Section III, where now the Raman excitation of 
one motional sideband is via the virtual excitation of the 
electronic transition with the lower state \g). The en- 
gineered Kerr operation V = e^ ** e -i<f>i**n e -icf>2cr x n 2 
generates all the circular coherent states in Eq.QJ 
of the vibrational degree of freedom for <j>2 = tt/M 
|35l | . In particular, we can see that the compass 

state (M = 4) is generated because e ±l ^ '' IV)n \a) = 
1 ( e ±*(V4) (| a j _|_ a )) + (| io ) + I _ The total 

unitary operation is U = V(ir/4)D c (a), and the out- 
put state is of the form given in Ea. (|24(l . with the 
electronic part as | |) = | |x) an d I i) = I ix)> and the 
vibrational motion of the ion in compass states of the 
form 



cat4> = 



cat 4 



™/ 4 (\a)-\-a)) + \ia) + \~ia) 
— a)) + \ia) + I — ice 
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Here a = ae 1 ^ 1 , and a = ae 



These two compass perturbed state is 



states are rotated with respect to each other in an angle 
2<fii. If it is desired to have the compass states equally 
oriented in phase space, as in the first approach, a con- 
ditional rotation could be performed. 

The engineering of the interaction Hamiltonian can 
be accomplished with N = 2 pairs of Raman lasers 
(with Lamb-Dickc parameters r/j and Rabi frequencies 
Qj, j = 1,2) driving resonantly the carrier (k = 0) 
transition. Indeed, one can choose the coefficient of the 
quadratic n 2 part in Ea. (|18|> to be A2 — 1, and impose 
that the cubic coefficient identically vanishes, A3 = 0. 
This sets the time t* when the compass state is gener- 
ated as <f>2 = \QL\A2t* /2 = 7r/4. The compass state 
is formed at the time t* irrespective of the values that 
the phases 0o and <f>i take. Thus, we do not require 
any fixed values for Aq and A\ in the interaction Hamil- 
tonian, simplifying the engineering design to only two 
pairs of Raman lasers. Selecting the set of Lamb-Dickc 
parameters r)j for the two pair of lasers as {0.4, 0.35} 
yields for the relative Rabi frequencies flj/flL the values 
{—520, 1154}. This requires that the Rabi frequencies of 
the two pair of lasers be related as | O2 1 ~ 2.2|f2i|, and 
have a it phase shift with respect to each other. The coef- 
ficients Aq and A\ are dependent quantities that can be 
obtained from Eas. l|16ll7|) . and are equal to A R* 605 
and A\ f=a —57. In order to minimize t* we choose 
the highest experimentally available Rabi frequencies, 
|Ox I ~ 5MHz and |0 2 | ~ 11MHz. Therefore, the ref- 
erence Rabi frequency H,l is ~ 9.5 x 10 3 s~ 1 , and 
the time t* at which the compass states is generated is 
equal to t* = 7r/2|f2i| w 165/is. The total time needed 
to generate the final state in Fig. 6 is equal to the 
duration of the conditional displacement (approximately 
equal to 10/iS for \a\ ~ 3, see previous sub section), plus 
the time t* for the compass state generation, that is a 
total time of approximately 175/is. 



V. WEAK FORCE DETECTION SCHEME 

Once the intermediate state of the form Eq. 1)24(1 is cre- 
ated, either by means of the first or second approaches, 
the probe is subjected to the perturbation to be mea- 
sured. We are interested in weak classical forces that 
couple to the center-of-mass motion of the ion, and thus 
cause a small linear displacement of the motional quan- 
tum state of the ion, irrespective of its internal electronic 
state. This can be simulated by applying the displace- 
ment operator D c ((3 = e lip sa/\a\) on the vibrational 
state of the ion with two pairs of Raman lasers, each 
of them inducing an equal blue sideband transition on 
the \g) and |e) states (see Section III). Alternatively, a 
uniform electric field oscillating near the trap frequency 
results in a displacement perturbation acting equally on 
both \g) and |e) [3(|. The general form of the created 
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|$) = -j= |cat 4 ( S )) I T) + ^= |cat 4 («)> I I) • (25) 

After reversing the evolution along the lines described 
in Sections III and IV we get the final entangled state of 
the form in Eq.©. In the first approach we have 



P g \g,n) = & \0)+A 2 \^a)+A* 2 



2 l")) ® 1.9) , 

(26) 

Ai(|a|s) = cos (b + (ip)\a\s)x 

A 2 (\a\s) EE ( e -i2sin(v0|a|o _ 

e i2cos( v )|a| S )/ 4) with b± ^ = CQS ^ igin^). Hence, 

neglecting contributions of the order C(e~' Q ' ) the prob- 
ability to measure the ion in the internal state \g) is given 
by 



where a = 2ae l 4 and A\ 
x cos (b-(<p)\a\s), A 2 = 



P g (s) = Al(\a\s) 1 + 2 



as 



(27) 



where we recognize in A 2 (|a|s) the approximation in 
Eq.Q for the fidelity function f(s) = \ (cat4|cat4(s)) | 2 . 
Thus, P g (s) exhibits the characteristic oscillatory behav- 
ior with a typical wavelength ~ l/|a| that allows, after 
inverting the relation in Eq. l|27|l . the measurement of s 
with Hcisenberg limit precision. 

In order to invert the function P g (s) to obtain the small 
displacement s we have to know the prior information 
that < s < sq, where sq, given in Eq.JSJ), is the first 
zero of Pg(s). This is not a restrictive condition since 
we have the flexibility of setting up the value of |a| in 
the experiment in order for so = 7r/2b(ip)\a\ to be an 
upper bound of the expected displacement. A good es- 
timation of the unknown parameter s requires repeating 
the measurement several times. The uncertainty in the 
determination of the true parameter s can be estimated 
if we observe that after R repetitions the probability that 
the outcome \g) is obtained r times is given by a bino- 
mial distribution, that in the large- R limit can be ap- 
proximated by the Gaussian distribution in the variable 
£ = rjR, which can be regarded as effectively continuous 
IIOL l37l . In this limit the probability distribution for the 
(£ = r/R) is 



estimator s = P g 1 



Pis) 



1 



V2ttAs 



S2 



where the uncertainty of s is 



As 



y/(l -_P 9 )P 9 
Rial 



dP n 



dy 



(28) 



(29) 



where P g = 



P g (s) and y = \a\s. Hence, we see that 
we reach the Heisenberg precision for displacement since 
\AR|a| ?s V Rn and Rn is the total number of photons 
used in the measurement. 
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In the second approach we have 

</P a \g,*%) = (B 1 \0) + B 2 \-2a) + B 3 \(i-l)a) 

+B A \-(i + l)a))®\g), (30) 

where 

B\ = (cos(aj"|a|s) + cos(aJ|a|s) + 

+ cos(a~|a|s) + cos(a~|a|s))/2 , 
B2 = (cos(aJ|a|s) — cos(aJ|a|s) + 

+ cos(a~|a|s) — cos(a~ |a|s))/2 , 
B 3 = i(e- 4( T+l Q l 2 )(sin(a+|a|s) + sin(a"|a|s)) + 

+e 4( *- |a|2) (sin(a s -|a|s) + sin(a+|a|s)))/2 , 
Bi = ^(e- ^( ^- |a|2) (sin(a+|a|s) -sin(a-|a|s)) + 

+e i(f + I°l a )(sin(a7|a|a) - sin(o+|o|a)))/2 , 

and we define af = 2sin(<£ ± <f>i), a~ = 2cos(ip ± 0i). 
In this case, the probability of measuring the ion in the 
internal state \g) is 

p g( s ) = \[ l + (cos(a+|a|s) +cos( a -\a\s) + 

+ cos(aJ\a\s)cos(a^\a\s))/2] , (31) 

that also oscillates with a typical wavelength ~ 1 /|o;| - 
Hence, following the steps described for the first approach 
we see that after inverting the function P g (s) we get also 
in this case the displacement s with Heisenberg limit pre- 
cision. 



the vibrational degree of freedom of a trapped ion are of 
the order of 100ms (see ref. [2^ an d references therein), 
so that the typical decoherence time for the vibrational 
compass state with \a\ ks 3 would be od the order of 
10ms. This decoherence time-scale should be much larger 
than the total interaction time for weak force detection 
(generation of the compass state, application of the per- 
turbation, and inversion of the dynamics), whether the 
first or the second approaches for compass state gener- 
ation is used. Assuming that the duration of the dis- 
placement perturbation takes approximately 3/iS (com- 
patible with the perturbations used in j3|| for engineer- 
ing the ion's reservoir), the total interaction time using 
the first approach is around 150/is, and using the second 
approach 353/is. We see that the typical decoherence 
times are much larger than the total interaction times 
in both approaches, pointing towards the experimental 
viability of the proposed scheme for quantum-enhanced 
measurements using the motional state of a trapped ion. 

The total interaction time for the compass state (M = 
4) generation is shorter when the first approach is im- 
plemented. It is worth noting, however, that the same 
tailored Hamiltonian implemented with the second ap- 
proach generated higher order (M > 4) circular coherent 
states for shorter times, as the necessary generation time 
scales as <j>2 = ir/M. Such higher order superpositions of 
coherent states also have sub-Planck phase-space struc- 
tures, and in principle could also be used for Heisenberg- 
limitcd quantum metrology. 



VI. DISCUSSION 



The scheme proposed in this paper for Heisenbcrg- 
limited sensitivity to perturbations with continuous vari- 
ables relies on the creation and manipulation of meso- 
scopic superpositions states of the motional degree of 
freedom of a trapped ion, that possesses sub-Planck 
phase-space structures. Any decoherence process, such as 
amplitude or phase decoherence affecting such quantum 
superpositions [Sfil ]. may destroy such small scale struc- 
tures, limiting the usefulness of the method for quan- 
tum metrology applications. Typical damping times of 
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